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ABSTRACT 


ihe “objiective of the project described in this report 
was to develop computer systems which would generate the 
euement Connectivity, and modal point co-ordinates for two 
and three-dimensional finite element programs using isopara- 
metric finite elements. The computer systems and sample 


\ 
problems are discussed. 
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i.) SENTRODUCTION 


The increased popularity of the Finite Element analyti- 
Cal method due to the availability of high speed, large 
memory computers has led to the solution of many heretofore 
unsolvable problems in structural analysis. Although users 
Se chas method welcome the greater freedom given’ them to 
pursue more fundamental problems, they are presented the 
tedious task of generating and checking the mesh input data 
Meeenewr particular computer system. “This data includes the 
arrangement of nodes within each finite element (connecti- 
Wity) and the co-ordinates of all the nodes. 

Automatic mesh generation is an attempt at simplifying 
input data for the finite element programs. If the user 
can communicate with the computer in a better, faster, and 
more efficient manner, money and labor will be substantially 
Sayed, and the possibility of human error will be greatly 
reduced. | 

Many automatic mesh generation schemes exist, however 
these srrpilenie are normally proprietary, very specialized, 
and not generally adaptable to finite elements having side 
nodes. The present system is an attempt to remedy that 
Srevation. 

Reterences, 3 and 6 prompted this "author to devise a code 
for two and three-dimensional automatic mesh generation 


Schemes, Utilizing the technique of mapping from a. local 


non-dimensional system of reference into the global cartesian 
Eyocen Of reterence. 

Reference 4 is a previous work on a three-dimensional 
automatic mesh generation scheme using isoparametric elements, 
however the coding incorporated an interpolation scheme which 
ted) to ili-conditioned elements in Some instances. 

To ailustrate the effectotereduction in the amount of 
@ata as a result of automatic mesh generation, the following 
example is introduced: Figure l(a) is a rather complex 
geometrical shape in the x-y plane, which can be sub-divided 
meee sections (super-elements) forswhich one of the i1sopara- 
Metric elements is avery close geometric approximation. 
These super-elements may be further discretized into smaller 
elements. The connection of the super-elements and elements 
iemore casily illustrated in a local non-dimensional system 
of co-ordinates, as shown in Figure 1(b), where all the super- 
Ewements are Squares. With the appropriate boundary descrip- 
tween Of the super elements as shown®in ’Pigure 1(a), and a 
super-element connectivity diagram as shown in Figure 1(b), 
enue fine grid of Cipune l(b) can be mapped back into the 
original object GO Produce ene w¢esired Mesh OL Figure 2. 

Four automatic mesh generation schemes will be presented 
mi this thesis. One is for two-dimensional linear, quadratic 
and cubic elements. The others for three-dimensional linear, 


quadratic and cubic elements respectively. 


Ga) 


Connectivity of super-elements in, Tocal system 
1(b) 


Fagure I. SEL. Dideramsr. 


Puplre (2. Final Mesh’, 
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Il.  ITSOPARAMETRICSEINIULE ELEMENTS 


The 41Soparametri¢e concept willbe discussed in.this 
ehapter in suffLicientidetail, to,establish only the procedures 
fom CONStructing isoparametric elements and the validity of 
ZetatirOnshiips used*to obtain a°*mappang from the local non- 
dimensional reference into the global cartesian system of 
reference. The readers are referred to references 1 and 2 
fore turther details of the isoparametric concept. 

The serendipity family of isoparametric elements [1] was 
mealuzed im this paper due to the ease of transformation 
(mapping) from a local non-dimensional system of coordinates, 
where all elements are squares or cubes, to a cartesian 


system with curvilinear boundaries. 


A. TWO-DIMENSIONAL ELEMENTS 
Consider the square isoparametric elements of Figure 3. 


Figure 3(a) is a linear element consisting of four corner nodes. 


7 Bie .\6 
linear element quadrac1ce element cubile element 
3(a) 3 (b) 3(c) 
Figure 3. Isoparametric Element Configurations 


Hie Locale.) Coordinates. 


| 
] 
y 


* 
. 


Figure 3(b)) 2S,a quadratic ellement consisting, of fourscoruner 
nodes and four side nodes. Figure 3(c) 1s a cubic element 
consisting of four cormer nodes and eight side nodes. 

A. curvilinear element, -sucheas that: ofo Figure 4,.may be 
produced by mapping the appropriate square of Figure) 5 into 


the x-y plane using the xeliataons 
n n 
ae ae (1) 


where (x53) are the cartesian coordinates of node i. 
N, are/Shape functions of the Serendipity family, listed in 
Section II.C.1., n is the number of nodes on the element 
finer ,/4 for linearselements, .8 for quadratic elements, 12 
for cubic elements). 

Note that the convention for node numbering within any 
element, (see,Figures*3.and 4) .is to begin at (&3n) =) (b,)) 


and to number the nodes proceeding counterclockwise. 


Figure 4. Curvilinear Element Mapped 
imsGartestan Coordinates. 
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B. THREE-DIMENSIONAL ELEMENTS 

Linear elements, Figure 5, have eight corner nodes. 
Ouadratic elements, Figure 6, have eight corner modes and 
twelve side nodes. Cubic elements, Figure 7, have eight 
Eemner modes and twenty-four side nodes. 

Similarly, as with two-dimensional elements, three- 
dimensional curvilinear elements may be produced (mapped) 


moane the relations 
\ 


where (X55 22;) arey_the cartesian coordinates of node 1. 

mes Shown in Figures 5, 6 and 7, the convention for node 
wuepering within any element 15 @to begin at (é,n,%5)— (1,1,1) 
and to number the nodes proceeding counterclockwise about 


tae zeta axis. 


Ooi a1) 


#5 (jas 1) 


Sol, -15,-)) 


Figure 5. Three-dimensional Linear Element. 


Figure 6. Three-dimensional Quadratic Element. 


Figure 7. Three-dimensional Cubic Element. 


ie 


Gee SHAPE FUNGELONS 


eG as = €€. ERP MesnniiiMim, co See Semraiee 


1. Two-Dimensional Elements 
a. binear elements 


See Figure S(a) for identitication of the nodes. 


N; = (1/4) (E65 +1) (ny +2) (3) 


where EG = +] and Na sh Oe 


b. Quadratic elements 


2 


[See- Faoure= sb) for identifacation of the nodes. 


forner nodes: 1°,35.,557 


N, = (1/4) (1#€,) (+n,) (E,*n,-1) (4) 
where a = +]1-and as ca | 


Bade nodes: .2,4,6,8 
Nie tL) 2 = 2 ikea (5) 
for nodes 2 and 6, where oma zal 


N, = (1/2) +8,)-n*), (6) 


for nodes 4 and 8, where Ef 8 ills 


¢. Cubic elements 
see Ficure Sc) fom adentitication of nodes. 


Comer nodes: 1,4,7,10 


N; = (1/32) (1 +65) +n, )9(E7#n7) - 107 (7) 


rs 


eadernodes: °5,6,11 512 
z LZ 
Nip D7 32) Glee (tm) Citar), 


where b5 aegis noe WAS 6 


mide nodes: 2,3,8,9 


N, = (9/32) (1tn,) (1-€7) (1+98,), 


where Es =i 3, ees Zl 


os Three-dimensional Elements 


a. Linear elements 


Seewikreuce 5 for identamweatron of nodes. 


Wea /e) (ie) ania ), 


where E3 = +1, Mea fans te) = A, 


b. Quadratic elements 


Scenrucure 6 for identiticacvon of nodes. 


Pomenenmodes: 1,3,5,7,13,15,17,19 


(8) 


(9) 


(10) 


Ny 1/8) (ee Cae Ee eng c 2) 
where eG =i eas sa) Us a pal 
Bide nodes: 2,6,14,18 

. , . 

N= hil seal are.) (12) 
where 25 = 0°, ee =a aes 1. 
Sige nodes: 4,8,16,20 

N; = (1/4) (l-n*) (140,) (+60) (13) 


Ale 
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GY» 


where a eal Ia ee Og Ae a= 
Bide. nodes: '9).10), 12 


N, = (1/4) (1-07) (1+€,) (tng) (14) 


where a Se aes ps i) eee coal ree ae) 


GaeCubie elements 


See Fapure 7 for adentativeation of nodes ~ 

Romer nodes: 1,4,7,10,21,24,27,30 
747°) -19] 
(15) 


N; = (1/64) (1465) (Lng) (+04) (9 (67 


where b3 = tle s tabs =) tl. 


erdeenodes:  2,3,8,9,22,23,28,29 


N; = (9/64) (1-67) (1#9€,) lng) (1+z,) (16) 


where S05 SAAS en) eee etl) Cs ee all 


1 


Pageemodes: 5,6,11,12,25,26,31,32 


N, = (9/64) (1-n7) (1#9n,) #55) (+E) (27) 


where a = +], ne > #1} 3, cane li 


Bmdernodes: 13,14,15,16,17,18,19,20 


N; = (9/64) (1-c7) (#94) +E.) tng) (18) 


where See #1, ean ars Gas aol We 


15 


Til. BASIS OF AUTOMATIC MESH GENERATION SCHEME 


As mentioned an the iantroductson, xreterence 3 prompted 
Pie utilization of the following scheme in this thesis 
eesearch. 

THe essence of the scheme is the use of simple mapping 
functions to transform simple orthogonal connectivity dia- 
mramssinto real geometrical shapes. Considering the partic- 
Mlar two-dimensional case of a parabolic quadrilateral of 
muswces4> answhach) the x and y co-ordinates of eight nodes 
meceKnOWwn, welGan write equation (1) from Chapter II, where 
N is a shape function associated with each node and defined 
mn terms of a local non-dimensional coordinate system &€ and n 
wiich has values ranging from 1 to -1 on opposite sides. 

ipacal shape funetions aré given im If.€., and for this 
wecmole equations (4), (5), and (6) apply. “If the co-ordinates 
Steehesnodal points are known, then cartesian co-ordinates of 
any specified point (&,n) can be found by equation (1). 

If the whole ae in which themeéesh is to be generated 
could be described adequately by a quadrilateral of the shape 
in Figure 4, a mesh of any refinement could be automatically 
Bemerated inside it by specifying the coordinates of the 
eight boundary nodal points and the number of sub-divisions 
mothe € and n«directions. The above scheme is shown pic- 


wOrially in Figures 8(a) and 8(b). 
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Carcular 


Duadratic Super-Element Quadratic Super-Element 


ihocal coordinate cartesian coordinate 
system system 
8 (a) 8 (b) 


AC 


Super-Element Selection 


Site) 
Buper-Element Type 
A quadratic 
B linear 
CG linear 
D Eubie 
iE quadratic 
wy corresponds to corner nodes 
© corresponds to side nodes 
Bh Corresponds to side nodes for a single Super-Element 


whose adjacent Super-Element is of a different type 


Fugue: 8 


7 


If the whole region could not be adequately described by 
Figure 4, then we could subdivide the region into more than 
ene quadrilateral and apply the same principles to each as 
mpecified above. These region subdivisions will be referred 
to as Super-Elements hereafter (see Figure 8). 

The local Super-Element of Figure 8(a) can be quite 
erastically distorted, even up to a point of making two sides 
mre along the same line in the cartesian system such as in a 
mesiole; Or Semu-circle. Care must however be taken not to 
make any corner angles greater than 180° as a non-uniqueness 
may result [1]. In addition, side nodes should be equally 
spaced between corner nodes to insure equal spacing of sub- 
movasions in the & and n directions and thus the x and y 
werectcions in the mapped cartesian system. 

ane criterva above also apply to the linear and cubic 
isoparametric Super-Elements for two-dimensional meshes, and 
ne Pinear, quadratic and cubic isoparametric Super-Elements 
for three-dimensional meshes. 

ihe decision on whether to use. lanear, quadratic or 
cubic Super-Elements will depend entirely upon the region 
boundary to Be mapped. Linear Super-Elements will map 
Guadrilaterals with straight sides; quadratic Super-Elements 
W211 map quadrilaterals with straight or parabolic sides; 
Giuac Super-Elements will map quadrilaterals with straight, 
BeeaDOlic. or cublre sides. Figure 8(¢c)° is an illustration of 
super-Element selection based upon region (structure) boun- 


dary conditions. 
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The final generated mesh for a given structure may con- 
tain nodes on the boundary which do not exactly ‘coimearde 
mith the actual structure boundary, or the elements within 
mayenot be represented sufficiently accurately. “In this 
case it may be necessary to have some of the coordinates 
majusted once the mesh has been generated. This adjustment 
will be quite easily accomplished with the aid of computer 
=oes, which can be obtained as output to the mesh generat- 


me programs discussed in Chapter IV. 
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IV. DISCUSSION OF COMPUTER PROGRAMS 


The computer programs presented in this thesis were 
mereten to SUppOoTrt the computer Systems “PLISOP' and “TRISOP': 
both of which were coded by Professor G. Cantin, and his 
Seudents of the Naval Postgraduate School, Monterey, Cali- 
mania,  TRISOP performs a structural analysis of three- 
dimensional problems using isoparametric finite elements, 
maale "PLISOP" performs the same analysis of two-dimensional 
meoplems. The input data to each system is quite extensive 
math regard to element connectivity and nodal point coordi- 
nates. 
| The automatic mesh generating systems discussed in thas 
chapter compute the element connectivity and nodal point 
eoordinates required for the above two computer systems. 

Pour mesh generating programs will be discussed. The first 
Mememaces data for 'PLISOP' with linear, quadratic and cubic 
elements, and is called 'PLIMEG'. The other three programs 
generate data for “TRISOP! With linear, quadratic, and) cubic 
elements respectively. They are, called RIMEG—i'. )  TRIMEG=Z ' 


' 


ana ‘TRIMEG-3'. 


A. THE GENERAL PROGRAM 
while” a “simple mesh generator using only ome super-— 
element may be sufficient for some purposes, it 15 limited to 


Smiple geometries WLth but a single unitorm property... For 
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complete. generality, as.suggested at the end of Chapter III, 
ene scheme was extended to many super-ellements, exther linear, 
Muadraticl, Cubic,;or any combination of the three, and repre- 
mmcedin.a)-'checkerboard' pattern, where each of the super- 
elements may be sub-divided and each may define a material 
With a single property. The above scheme is presented in 
Sections B and C of this chapter for two and three-dimensional 


Meanes respectively. 


B. TWO-DIMENSIONAL MESH GENERATION (PLIMEG) 
Pee DeLinttrons “(see Figure -9) 


an oii eran acronym; used, #60, meter Co a suped- 
element as defined previously on page 7. 


ba Master Grid: A 10 x 10 “checkerboard” array 
Of SELS which is utilized as a working 
board for SEL numbering and orientation. 
The master grid’ SELS are numbered sequen- 
tigtly alone grid columns: 

ce Rows: ~SUbdavisions of a SEL in 1 darectaon 

de Columns: ) subdivisions of a SEL im & direction 


er Lklenments: those portions Of 4 Sklascreated by, 
rows and columns 


Ze ouper-clement Selection 

The number of SELs required in a given mesh is deter- 
mined from boundary considerations and the types of materials 
from which the object is composed. Enough SELs must be 
incorporated to adequately define the boundary of the object, 
and at least one SEL must be utilized for each different type 
of material. Figure 10 shows an example for a particular 
ebyect: The redder will note in Figure 10(a) “that the 


common boundary of SEL 1 and SEL 2 is based on material type, 


aM 


SEL f12°(8 nodal point) 
Master Grid 9(b) 


9 (a) 
Avindicates corner nodes; @ indicates side nodes 


All nodes indicated comprise SEL boundary nodes- 


Fugure 9  PLEMEG Definitions. 


While the remaining boundaries are determined from geometri- 
cal considerations only. 

Once pthe quantity and. typeror skus have been deter 
Mined, the proper SEL orientation in the local nondimensional 
meoeemeof rererence must be decided upon, as shown in Figure 
we (b)i. The local system is now superimposed upon the master 
grid of Figure 9 to obtain the SEL numbering scheme of 
Figure 10(b). Note: this superimposition may be anywhere 
Within the boundary of the master-grid. 

oe 6Mesh. Createria 
the following rules must be adhered to at.all times: 
a. All elements within a given mesh must be of the 


mie eypes(t.e., laineay, Or quadratic or cubase), ‘and only 


ene type of SEL May vary: 


2 fh 


oe 


Material 


Material 2 


1 


Cartesian SEL Diagram 


EO! Gay) 


Gorresponds to corner nodes 


corresponds to side nodes 


eorresponds to side nodes for 
ausungike SEL whose adjacent 


Spe 1s of a different type Be eee ee 


10 (b) 
m4 rows, 5 columns (linear) 
aes rows, 5 columns (quadratic) 
i> rows, / columns (quadratic) 
iS rows, 6 columns (cubic) 


Final Mesh 
MOE) 


Figure 0. PRLIMEG Example. 
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b. The number of rows in all SELS of any master grid 
mow must be equal (ive., SELS Z,12,22...92) 

¢. The number of columns) 1m all SEs of any piven 
Geter grid column must be equal (.e., Skis Ti 12,15. . 220). 

d. SELS may be Connected to adjacent SEES 1m the 
following ways (see Figure 11): totally connected, connected 
at one corner only, totally disconnected, and connected at 
menace corner nodes only but not along common edge. Figure 11 
euows only horizontally connected SELS, but the same connec- 
fion- schemes apply to vertically connected SELs as well. 

4. Bandwidth, Node and Element Numbering 

Element numbering is always along master grid 

columns, and with the exception of discontinuities between 


Mimaeecent oLLsS, such as in Figures li(b)” and 11(d), node 


nN & 
, AN 
Totally Partvaliy Totally 
Connected Connected Disconnected 
a Cayie: eG Gow) IGG) 


Connected at 
Corner Nodes Only 
VACd) 


Prune niin. okl, Connections. 
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numbering is also along master grid columns. Therefore, to 
obtain a node numbering scheme which results in the best 
eanputacaonal efficiency, since "PLISOR” and "“TRISOP “ware 
designed to be used with a banded stiffness matrix solution, 
meets prudent that the number of rows in the master grid row 
Sontaining the most columns be less than or equal to the 
total number of columns in that particular master grid row. 
mine reader is noeaned to Section Casson chs chapter for 
additional information regarding this topic.) For example, 
eae rFigure 10 the master grid row containing the most columns 
eocrd GFOW tWO, im which there are sixteen columns, while 
the number of rows in this grid row is five, thereby satis- 
Pomp ethe above criteria for bandwidth. Another example 
concerning bandwidth is as follows: suppose the object to 
memconsidered a5 in the shape of an “L" and that all SELs 
Seneain an equal number of rows and columns (i1.e., 3). Two 
eeorerent SEL combinations could be utilized to describe 
Ene mesh as shown in Figures 12(a) and 12(b). However it is 
eyadent from Figures 12(c) and 12(d) that bandwidth would be 
Peecater for the shape of Pipuie “bZiGbi,sameec: there ane Six 
rows in master grid column one for Figure 12(d) vice three for 
Puegure 12(c). 

There will be many instances when the geometries of 
feeure 12 will present itself, and when this condition arises, 
the user must determine whether element shape within the mesh 


Or a mesh with minimum bandwidth is the more crucial. 


Z5 


Cartesian SEL Diagram Cartesian SEL Diagram 
12 (a) U2 Gby) 


ee ee ee ee a ee) 


Local SEL Diagram 
ANS) 


Local SEL Diagram 
12 (d) 


Bogune 12. The “LY Shape. 


Whenever discontinuities such as those of Figures 
f(b) and 11(d) exist in a mesh, a good guideline to follow 
ms to have them exist above or below a SEL, if possible, 
Mmaeietm chan to the lett or right, since the latter type of 
sel connection occasionally results in a somewhat non- 
efficient node numbering scheme due to the complicated 
Podin: ImecOorporated in Subroutine CONN of (PRIMEG:. 

seo EP LIMEG Output 

The major portion of the output wath, PLIMEG as 
element connectivity and nodal coordinates, however addi- 
raonal output has been incorporated to aid in the use of 
feasOP.” To begin With; all data input to PLIMEG will be 
reproduced as output as a check to the user, in event a 


mispunched data card is experienced. The total number of 
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elements and joints in the mesh together with the half 
bandwidth willialso be printed. «In addition, an option to 
punch on data cards the element connectivity and coordinates 
Gr the joints 1S incorporated. Lastly, a plot of the mésh 
may be obtained on the printer, whereby Subroutine GRID calls 
eeoroutine UTPLOT in the IBM 360 source library. However, 
since plots on the printer are somewhat distorted and only 
SausiSt OL points, a plot can also be obtained on the off- 
Panes plotter (Calcomp) whereby Subroutine GRID calls ,Sub- 
feucine DRAW invthe source library.  The-resulting+» plot .is 
Bsethose Shown in Figures Z and 10(c). The printer plot 
should be used in determining proper mesh discretization and 
euen Calcomp for the final plot of the mesh. Both plots can 
not be obtained on the same computer run. 

Appendix A discusses input data preparation for 
PLIMEG while Appendix B discusses sample problems and shows 
the input/output formats, the node and element numbering 
mevemes, and the resulting plots. These two Appendices may 
be removed and utilized as an instruction manual for PLIMEG 
once the contents of this and the previous chapters are 
comprehended by Ehe suse. 
C. THREE-DIMENSIONAL MESH GENERATION (TRIMEG) 

PRIMEG actually consists of three separate programs: ‘one 
memeecach type of element (1.e., linear, quadratic, or cubic). 
imey are called TREIMEG-1, TRIMEG-2, and TRIMEG-3 respectively. 
However, since the computer scheme and the data input incor- 


porated in these programs is identical, they will be discussed 


7) 


as. one (i.e., TRIMEG). TRIMEG alsomimecorporates, many of 

tne ideas of PLIMEG, therefore only the .dissamrlaritzes 

between the two schemes will be presented in this section. 
Pe Vetini trons (see maeure 13) 

a. Master Grid: a 5x5x5 three-dimensional "checker- 
board array of SELs. The master grid SELs are numberéd 
sequentially along master grad Columns sim Succes save=master 
Brageslices. 

Db. Slices 3) qsub-divastons of -SELein® G4 direction 


Gc. Element: “those portions Of al SEL icreated by 
rows, columns and slices 


femeeiesh Criteria 

ihe following rules must’ be vadhered to at all times. 

a. the munber of rows ingall’ SELS of amy master grad 
PMnct De equal (1.e7, 256. .«°21,26,51....46551,56...71,760, 81, 
meee. O0l ,106,..121). 

db. The-number of  onTNs inea li SEES@on tanyimas ter 
maancolumn must be) equal (ive., 1-5,20-30,51-55,76-80, 
POT-105).. 

eG, ehesnumber of Slices ini alieiSELS tor any master 
Metceslice musit be equal (i.e., 26,2728: ...50): 

Ge olhertaees sot tadjacent SELsamayebe conmected in 
aay on the ways shown in Figpune ll (applies to all six faces), 
Peeme the exception lof Figure lil(d). If this situation exists 
am amesh, then four SELs must be utilized as shown in 


Faeure 14. 
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Figures 14;  TRIMEG SER, Connection for 
Samidar Connection of 
Freure Lian 


3. Bandwidth, Node and Element Numbering 

Node and element numbering 1s identical to that of 
PiiteG Lor each-successive slice in the mesh. The following 
rules apply t. obtain minimum bandwidth. The number of rows 
in the master grid row containing the most columns should be 
less than or equal to the.total number of columns in that 
articular master grid row. In addition, the number of 
columns in this master grid row should be less than or equal 
fomene’ total number of slices in this particular grid row. 
This scheme can be visualized by sketching a single SEL with 
one row, two columns ane three slices and then numbering the 
nodes along columns in successive slices. The half-band 
faaten 15 50 as shown in Figure 15(a). Using two rows, one 
eolumn and three. slices, produces a halt bandwadth of 33, 
eee ihe reader, should sketch several possibilities for 
few, COlumn and’silice combination (Figure 15) to convince 
himself that orientation of the mesh determines bandwidth, 
and thus efficiency, when solving problems with PLISOP and 


ERISOP. As can be seen in Figure 15, numbering is itdentical 
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half-bandwidth 
15(a) 15(b) 


12 


5X2 LXZxS 
half-bandwidth = 54 half-bandwidth = 36 
15(c) Ts (a) 
Breure. 215i: 


op 


eo that, of PLIMEG for 2 sinele SEM wathwonecmolace: wietnour 
the third dimension. For PLIMEG the half bandwidth for 
Fipure 15(a) would be 8. 
4. TRIMEG Output 

TRIMEG output, is, rdentical ito) PLIMEGeexcept formmtme 
piottane package. -With TRIMEG a plot can only be obtained 
on the Calcomp plotter. The plots are of three different 
prowections on the planes produced by rotating the mesh 
about ehe: Z-aX1s, the fattest x-axis and the ilacest y-axcus 
meapectavely. The. three angles of rotation (Euler angles) 
and the resulting axes are shown in Figure 16 and the scheme 


Pomcontained in subrowtine TRFR of TRIMEG. 


i] 


anele of srotation. about “z-axis 


angle of rotation about X,7axis 


angle of rotation about Y7axis 


Live monet xeS ROtAtLONS for Projection Plotting. 
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Appendix C discusses anput data preparation for 
TRIMEG while Appendix D discusses sample problems and shows 
the input/output formats, the node and element numbering 
schemes, and the resulting plots. These two Appendices may 
be removed and utilized as an instruction manual for TRIMEG 
once the contents of this and the previous chapters are com- 


prehended by the user. 


oY) 


V. CONCLUSIONS 


ihe wtrli zation of PLIMEG and TRIMEG automatic mesh 
generation schemes along with existing finite element pro- 
grams, now makes it much more practical to solve realistic, 
couplacated problems, since the amount of input data has 
been reduced as much as 95 per cent. This decrease in the 
amount Of input: data’ cards required for the finite element 
programs not only reduces the chance for error, but makes 


problem solving more of a challenge and not just a headache. 
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VI. RECOMMENDATIONS 


PEGT OF STRUCTURE MESH 

In addition to the present plotting routine incorporated 
in 1RIMEG, it would be desirable to obtain a plot of the 
generated mesh which did not contain hidden lines. Profes- 
sor Gilfes Cantin of the Naval Postgraduate School, Monterey, 
pees a is presently working on such a solution. Adding 


tne actual node numbers to the computer produced plots would 


also be Of great value. 
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APPENDIX A 


INPUT DATA PREPARATION FOR 'PLIMEG' 
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APPENDIX B 


PLIMEG EXAMPLES, INPUT AND OUTPUT 


Thais appendix discusses two example problems. The first 
will show SEL selection and connection procedures dealing 
Miecihean eight SEL object with discontinuities along some 
edges, and also the node and element numbering schemes incor- 
porated in PLIMEG. The second will be a more straightforward 
example illustratine the detail of input data and output 


Formats. 


Pee xAMPLEW“ONE™ (Padeye; see Figures 17 and. 138) 

ihe common!’ boundary of SEL Zifand SEL 22 a's (based on 
Material type, while the remaining SEL boundaries are deter- 
mined from geometrical considerations only. As noted in 
Figure 17(a), the common, boundary of Shi Zand SEL 42 musé 
Bewconnected manually once the mesh has been generated. 
Titsess accomplished as folliows' (see Figure 17(c)). The 
mmatemesh of, Figure 18 will not contain nodes 77-80, there- 
fore remove these nodes from the output punched deck of 
coordinates and change the element connectivity cards for 
emenents? 55,50;and 57. In other words, node! 77 (becomes 
node 1, node 78 becomes node 2, node 79 becomes node 3 and 
node 80 becomes node 4. The number of nodal points (joints) 
will now be 76 vice 80, while the number of elements will 
memain the same (a.e., 156). Halt bandwidth will be increased 
from Z0ece Las, but may ber held te 50 af SEES 2,254 5 lcales 


and I4, for example, are used to describe the mesh. 
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Material 1 


Connected 


7 (2) 


ms fos seman 
12 16 18 
Local SEL Diagram Showing SEL Numbering and Corner Nodes 
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Subdivided local SEL Diagram Showing Node and Element 
Numbering (linear element) 


7 Ce) 


Figure 7. 
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Figure 18. . Final Mesh for Example One. 


ihe remaining variables for input to PLIMEG will be shown 


in the more straightforward example which follows. 


B. EXAMPLE TWO (See Figures 19, 20 and 21) 

Figure 19 shows a simplified mesh generation scheme con- 
Sisting of two SELS, six elements (linear) and 12 nodal 
points (joints). Although SELs connected vertically could 
have been utilized, this would have established a mesh with 
three rows and two columns resulting in a half-bandwidth 
greater than that of the existing mesh. 

Although the mesh is quite coarse, it adequately describes 
fie data input for the given boundary situation. The SELs 
Eoula have been subdivided to produce a much finer mesh, 
however the amount of input data would remain unchanged. 

Figure 20 shows in detail the input and input format 
tequired for PLIMEG as described in Appendix A, while Figures 


Panda 22 present the output and “output format. The last 
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two columns of the connectivity matrix (Figure 22) are KIND 
jac TYPE. These quantities represent the kind of element 
anaeene type of material for that element, respectively. 
KIND equal to one implies linear elements, two implies 
Guacnacic, and three implies cubic. Whis information as 


Tequired input for PLISOP. 
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Figune 19 -Exampile Twe-. 
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pusures; 20. PLIMEG Input Data (example two). 
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PLIMEG Output (input data check). 
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APPENDIX C 


INPUT DATA PREPARATION FOR 'TRIMEG' 
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APPENDIX D 


TRIMEG EXAMPLES, INPUT AND) OUTPUT 


Since the covicept for TRIMEG is very similar to that of 
PRIMEG, only the figures illustrating two examples will be 
Showm in this appendix. The user should be completely 
familiar with the contents of Appendix B (PLIMEG Examples, 
inpuceand Output) prior to utilizing this appendix and 


Appendix C (TRIMEG Input Data Preparation). 


A. EXAMPLE ONE (Crack Propagation Problem) 

Figure 23 shows the local SEL diagram indicating SEL 
numbers and SEL corner node iiumnie ae The cartesian SEU 
diagram is not shown since it would be identical to that of 
EncesoOcal diagram since the object is a cube. Figure 24 
shows the final mesh, where each SEL contains 2 rows, 2 
columns and 2 slices. 

Example two will completely illustrate all input, output 


and respective diagrams. 


-B. EXAMPLE TWO (Folded Plate) 

Figures 25, 26 and 27 show a simplified mesh generation 
scheme consisting of three SELs, 27 elements (linear) and 
emehty nodal points (joints). Although the mesh is quite 
Gearse, it adequately describes the data input for the given 


boundary situation. 
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Figure 28 shows in detail the input and input format 
required for TRIMEG as described in Appendix C, while 
Papures 29 and 30 present the output and output format. 

Each SEL contains one row, three columns and three 


s1ces. 


od 


10 


Local SEL Diagram 


Pirpure 22. 6 Crack Propagation Problem. 


SZ 
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LWearowne 
N Mi x Se 


rN Y 
ON | 


4 


Carcesian SEL Dirapram 


Pigure 25.7 Folded Plate. 
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Figure 26." Local’ SEL Diagrame 


5) 


Final Mesh 


Regure, 27. Folded Plate. 
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TRIMEG Output (Input Data Check). 


Figure 29. 
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COMPUTER LISTING (PLIMEG) 
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